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Abstract 


This paper describes studies of electrons between 25eV and 9.9keV 
in the solar wind. The transport of energy in the re*c frame of the 
plasma is evaluated and shown to be parallel to the Interplanetary 
magnetic field. The presence of electrons from solar events causes this 
energy flux density to exceed the heat flow due to thermal 
electrons. In one such event, the observations are shown to be consistent 
with the solar electron observations made at higher energies. When 
observations are made at a point connected to the crrth's bow shock by 
an interplanetary field line, a comparatively large energy flux along the 
field towards the sun is observed, but remains outwardly directed 

during this time interval. In either case H is found to be consistent 

II 

with measurements made on Vela satellites by a different method. These 

-2 2 

values, less than 1x10 ergs/cm / sec, are sufficiently low to require 
modifications to the Spitzer-Harm conductivity formula for use in solar 


wind theories 
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Introduction 

In this paper, we describe electron observations in the solar wind 
which cover the energy range 25eV to 9.9keV, that is % £rom energies 
characteristic of the plasma electrons to energies just below those 
normally observed in solar electron events. We show that transient effects 
can be important in determining the electron energy flux in the solar wind. 
We define the energy flux density along the magnetic field in the rest 
frame of the plasma to be E^, where the subscript indicates that -the energy 
flows along the magnetic field. We contrast this with H , the heat 
transport along the field which we assume to characterise the plasma in 
the absence of transient effects at higher energies. Using the constancy 
of the electron temperature observed both by Montgomery et al. (1968) and 
by the present experiment as justification, time averages of the electron 
observations over a period of 2 to 4 hours are used to construct the 
distribution function over the electron speed range 2960 km sec"* to 58,900 
km sec *, corresponding to the energy range 25eV to 9.9keV. Computation 
of the energy flux density E^ as a function of the upper limit of the 
velocities shorts that E^ in general does not converge to a unique value. 

The reason for the divergence of E^ is that the thermal electrons 
characteristic of the plasma, emitted continuously by the sun as a whole 
and making up H^, have characteristics which are not easily distinguishable 
from non- thermal electrons emitted from time to time by different means, 
since their spectra merge into one another. In addition to the transient 
effects of solar electrons (Lin and Anderson, 1967, Lin, 1970), which can 
make large contributions to E^ between ~340eV and 9.9keV> a net energy flux 
towards the sun is observed when. the detector is situated on a magnetic field 
line intersecting the earth's bow shock surface. 
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As a result of the influence of these transient effects at the higher 
energies, the shape of histograms of E^(E u ), where E^ is a fixed upper energy 
limit, might only partially represent the heat flux variability. The present 
observations indicate that for E <0.340 keV, when E (E , the heat flux has 
values consistent with those of Montgomery (1971a), and therefore support 
his conclusions regarding the plasma conductivity. 

Montgomery et al. (1968) have calculated the third moments of the 
electron velocity distribution function in the solar wind using detailed 
measurements up to 1600eV from Vela satellites. Determining the energy flux 
density in the rest frame of the plasma, E , they interpreted this at 
quiet times as the heat flux H^. The contribution of protons to the heat 
flux is about one thousandth of that carried by the electrons. 

More recently Montgomery et. al. (1971a, b) has reported a most probable value 

of the heat flux <H > ~ 0. 6-. 7x10 ergs cm sec . Such a low value •'s not 

II - 

consistent with the predictions of current solar wind theory, using the 
Spitzer-Harm conductivity formula (Montgomery et. al.,. 1971a). 

The measurements refer to higher bulk speeds (~300 - 325 km sec *) than 
those predicted by the theories (''*>£60 km sec *). These higher bulk speeds, 
presumably requiring the deposition of energy in the heliocentric range 
2-25 solar radii, should accompany higher heat fluxes than those presently 
predicted, which are already much larger than the experimental results, 
Hundhausen (1970) . One-fluid theories with modified conductivity which 
can predict close to the observed values, using the observational data 
as a boundary condition, have been reported (Cuperman and Harten, 1970, 
Whang, 1971, Wolff, Brandt and Southwick, 1971). Because of the connection 
with the transport parameters of the plasma, measurements of E^ are of 
interest, especially if carried out by different methods. 



Experimental 
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The observations described in this paper were made using the 
triaxial electron spectrometer on 0G0-5, described elsewhere (Ogilvie, 
et al. 1971). Three cylindrical electrostatic analysers accepted 
electrons whose directions of motion fell into three cones of ten degree 
half angle. The axes of these cones were mutually perpendicular , and 
made equal angles to the normal n to the mounting face of the satellite 
body. The attitude control system maintained the direction of n pointing 
always away from the earth. The data to be discussed here were taken 
when the spacecraft was in the interplanetary medium close to apogee, 

0-24 R g ) on outward legs of its orbit. The angle between the helio- 
centric direction and that of n varied less than one degree per hour 
during these periods, so that the direction of the detectors can be 
considered to have been fixed relative to the earth. Every 23 seconds 
the analyzers simultaneously measured the energy spectrum of electrons 
entering at three different angles {3 with respect to the magnetic field. 
The constancy of the detector sensitivity was checked by 3 radioactive 
sources and intercalibration between the detectors was checked using the 
data itself. 

Technique 

The basic assumptions of our analysis are a), the axial symmetry, 
and b). the near time- independence, of the velocity distribution about 
the interplanetary magnetic field line. Thus in a polar coordinate 
system, with the magnetic field direction along the polar axis, the 
distribution function is independent of the azimuth cp and the time t. An 
argument in favor of such an assumption is that T^,, a statistic of this 

distribution, is found experimentally to be essentially constant in time 
and the electron anisotropy is small (~1.1 - 1.2) even when proton 
temperatures are highly variable, (Montgomery 
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et al 1968). We therefore assume that <f (v,@(t) ,a(t))> constructs 
<f(v)>t, where v is the magnitude of the total velocity v; |3(t)*cos *(B(t)«V) 
is the angle that makes with the field line, <y( t) the angle between u and 
v as a function of time. In other words the time dependence of the measure- 
ments enter essentially through that of 8* since <y is almost a constant 
during a given measurement. 

Correlation Between Electron En e rgy Flux Direction and the Magnetic Field 
Direction in the Interplanetary Medium 

The energy flux vector of the plasma is the third moment of the 

distribution function, whose components 

Q j = JJf & || v-U|| 2 (Vj “Uj ) f (v , 8 , cp, t ) v 2 dvs ingdadcp 

* 2tt JJ* ^ w^wjf * (w,8* » cp* , t) Jw^dwsinB'dB 1 

where J s d(v,B,cp)/d(w,8* ,cp* ) is the Jacobian of the transformation into 
the rest frame of the plasma, in which w,3', and cp* are the coordinates. 

In these equations w * ^ - u, where v is the total velocity vector of the 
electron, u is the bulk speed, and w^ is the jth component of w . Thus, 
w is the random velocity with respect to the rest frame of the plasma. 

We also define the magnitude of the energy flux vector per unit solid angle 
as % 

0Q<a >Eu) = 


f* (w,$ , ,t)Jw 3 dw 


Where is the energy corresponding to the lowest velocity. 
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■ jE(w){f' (w,8* »t)w 8 /m}jdE(w) 

~ J*E(w){f 1 (w,0* , t) ^[1-21100 sa/v ]^/m}dE(w) 

where J :g ||W||/||Vi|. Since as f(v,0,t) 

f ' (w, @* , t)v 3 /m ss the measured differential flux. 

We thus obtain that 

Qq a Je(w)F1u*[ 1 - 2ucosa/v +u a /v*]3fedE(w) 

The bulk speeds used here were obtained from the MIT experiment on Explorer 
35 (Binsack, private communication), since they were not obtained from the 
present experiment with sufficient accuracy* The units in which this flux 
is expressed is particles/cmP/str/kev. It will be noted that Qq is a 
function of the highest energy employed in the distribution determination, 

Ey which is the upper limit of all the above integrals; the lower limit, E^ , 
is, unless noted otherwise, equal to 25eV. These and all subsequent integrals 
were evaluated numerically by interpolating a cubical spline 
through the data and integrating the result. Such an interpolant injects 
minimal prejudice to the quadrature determination in the sense that it is 
the smoothest curve which Interpolates all the data (Thompson, 1970) , (Scudder, 
1971). 

Let us examine <Qq> as a function of g, Three alternative possibili- 
ties for the energy transport are shown schematically in Figure 1. If there 
were no preferred direction of the energy transport with respect to the field, 
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<Q^0)> vs P would appear as shown In the upper panel. If there were a 
loss cone about the field direction such a plot would resemble the middle 
panel. The transport of energy essentially parallel to the field without 
loss cone (0»O or tt) would lead to a dependence of < Qq > upon P as shown in 
the lower panel, (or its reflection in (3 ■ tt/2) , where the flux is constant 
within a cone about the direction of B, (or -B). 

For several periods when the spacecraft was in the Interplane tary medium 
we have plotted <Q^(£)> in this manner. To avoid possible Interference from 
terrestrial effects, we have picked for our first examples times when the 
Interplanetary field had a low probability of bow shock intersection, and 
two such summaries are shown in Figure 2. On the right is an example taken 
from a positive magnetic sector when the field was directed away from the 
sun and U-380 km sec’ 1 , ard on the left from a "toward" or negative sector, when 
U-410 bn sec’ 1 . The data points are averages over five degree intervals of p. 
Whenever the p^ of. the ith detector of the triad is in the Jth angular interval, 
its value of QqOjCPj) * 8 averaged with other such determinations at subsequent 
times regardless of i. This technique directly exploits the assumed time 
independence and axial symmetry of f(w) as outlined above. Because of the 
variability of the field direction, more than one detector generally 
contributes to the average characteristic of an angular interval. All three 
detectors sometimes contribute to the high flux regions shown in Figure 2. 

If the forward cone is uniformly populated, as is indicated by a fiat step on 
the right, the interdetector calibrations can be checked, and small corrections 
to the data made by normalizing with respect to one detector. 

The March 8 example is from a negative magnetic sector, i.e. B field 
toward the sun. The upper (diamond) curve shows energy flowing away from 
the sun occupying a band of directions at about 120° from the magnetic field 
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direction, fpr an electron energy range of 25eV to 9.9keV. The lower curve 
(crosses) shows a similar graph covering the energy range 25eV to 340eV, 
When modified by the cos $ sin £ geometrical factor (see Eq. 3) this 
yields a small net heat conduction anti-parallel to the field, i.e. away 
from the sun. In the March 13 example two similar curves are presented 



the magnetic field direction was reversed, (note the change in vertical 
scale). The peak for the higher energies now occurs for angles <50°, and 
there is no empty cone. At lower energies a small energy flux awry from 
the sun is still measured. 

These plots are thus consistent with energy transport at all observed 
energies being parallel to the field line and flowing away from the sun 
regardless of the aveiage polarity of the sector structure. The minima 
of <Q^($)> in the transverse direction to the field suggest the inhibition 
of electron energy conduction in this direction. Figure 2 shows that 
the time average direction of maximum Q is correlated with that of B @ 
where £ is the earthward direction along an interplanetary field line. 

In Figure 3 we present similar evidence on a shorter time scale. Here 
Q^O(t)) is plotted isometrically against (3 and time for each detector, 
each point being separated by a time interval of 23 secs. These 
examples are from outward sectors and therefore when (3 is 
near 0, QqO) should be larger than when & is >70° • We have divided 
the energy range to E u into two parts; E^ and E u are 25eV and 340e? 

in the upper panel and 550eV and 2300eV in the lower panel. Both panels 
show that the detector having the smallest value of 0 measures the largest 
value of Q(P). Thus in each case the net energy flux E^ is closely 
correlated with the field direction ft. 
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Similar studies are being made of the transport of energy by electrons 

along field lines which intersect the earth's bow shock. Figure 4 shows an 

example obtained during a four hour period on March 19, 1968 when OGO-E was 

situated at X__*13.1 R -12.6 R , 2 13.7 R , very close to the bow 

&£ e sl e bL e 

shock in the dawn meridian when the interplanetary field sector was positive. 
Here <Q^> plotted against (3, and during this time field line intersection 
with the bow shock was almost certain. The crosses in Figure 4 represent 

calculated up to 210eV, and indicate a net heat flux H from the solar 

II 

direction. The diamond- shaped points, however, represent <Q^> calculated 
up to 9.9keV and indicate a large energy flux from the bow shock direction 
towards the sun. Protons of energy ~5 keV have been observed upstream of 
the bow shock (Asbridge et al. (1968), Frank and Shope (1968) which suggest 
that the present phenomenon is probably associated with electron flow from 
the shock as has been inferred by Fredericks, Scarf and Frank, (1971). In 
this example, all the observed backward cone angles contribute to the reverse 

flow, in contrast to Figure 2, where electron flux from the sun is confined 

o 

to a cone of half angle ~70 . Thus electrons from the sun were somewhat 

better collimated than those from the bow shock. 

The Magnitude of Energy Flow in the Solar Wind 

The energy transport parallel to the magnetic field is 

(E^) s Q.B « 2tt J n <Q n O,E u )> sii^ cosf3 d£ 

where the explicit dependence of E^ on E u is noted and the axial symmetry 

used. If transient energy flow from the sun as well as quasi-steady heat 

conduction, H , existed in the solar wind, E vs log E might look as in 
!! II ^ 
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Figure 5* In the upper panel the convergence of E below some energy E 

II c 

allows an unambiguous determination of H^»E (E c >. If transient energy 
flow occurs at lower energies, as in the lower two panels, H no 
longer equals E^(E c ) but is E^(E Q ) where E Q <E c . The bottom panel illustrates 
the case where the electrons providing the transient energy flow are 
associated with the earth's bow shock. The thermal electron population 
still provides a flux of heat away from the sun, and the population at 
higher energies a flux of energy directed towards the sun. 

Figure 6 shows E^ vs log E u plots for three periods of observations in 
the interplanetary medium, for which bow shock intersection of the threading 
magnetic field line is improbable. Each curve is from a different apogee 
on a succession of apogees in the interval March 8-13, 1968 and they are 
marked A, B, C on the inset diagram. The examples are successively more 
divergent, in the sense of Figure 6, as the day number increases. It will 
be seen that extrapolating the trend of the high energy parts of the 
curves to E^*0, indicates that the energy separating the electron heat flux 

H and the contributions to E due to transient electrons is about 340eV, 

II II 

since above this energy increases rapidly with E. Also, we see that at 

high energies E^ increases with time, whereas the value of (^E^(340eV)) 

remains about the same within the uncertainties of our determinations, 

-2 - 2-1 

and is certainly less than 10 ergs cm sec 

Comparison between the present results and data of Lin and Anderson 


(private comnunication) show that the highly anisotropic flux of electrons 
having energies between 340eV and 9.9keV tracks that of electrons having 
energies greater than 22keV. The electrons observed here are thus 
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associated with the same solar event as those observed by Lin and Anderson 
and are the cause of the divergences. This conclusion is made more 
plausible by Figure 7, where we see points from Lin and Anderson at 26keV 
and from the present experiment above 890eV obtained on March 8 and 11. 

The present observations are joined to those of Lin by the smoothest 
curve, (Thompson, 1970). At 26keV, the slope of the curve is that observed 
over the energy range 22-45keV by Lin. Although Explorer 35 passed into 
the tail before our next series of interplanetary observations, the Forbush 
decrease after the SC continued well past the March 13th observation 

t 


by OGO-E, which indicated further enhanced levels of fluxes in the non- 
thermal portion of the plasma distribution function causing the more 
divergent nature of E^. 

Figure 8 shows the value of E^ plotted against log E^ for a bow-shock 
Intersecting case. At energies below lOeV the energy flow, the plasma 
heat flux, is toward the earth, but at higjher energies the flow is much 
larger and towards the sun, away from the bow shock. The values of 
H =e^( 1600) presented by Montgomery were obtained by making detailed 
measurements of f(v) below E^ (+1600eV) on Vela spacecraft which were in 
nearly circular orbits of radius 19 R e » These orbits were closer to the 
average bow shock position than the regions where the present data was 
obtained. Since the Vela spacecraft do not carry magnetometers, the 
magnetic field direction at the spacecraft was not known, and partitionirg 
the data into sets associated with intersecting and non-intersecting field 
lines was approximately carried out by noting whether or not there was 
a significant component of heat flux coming from the earthward direction 
along the field. In addition, if outward streaming protons were present 
it was assumed that the spacecraft was on a bow shock connected field line, 
(Montgomery - private communication) . 
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In view of the above discussion of transient effects due both to bow 

shock Intersection and to solar electron production, Is not always easy 

to deduce from such measurements. It Is therefore very satisfactory that 

-3 -2 -1 

the values of H^(~4xl0 ergs cm sec ) obtained here are In good 
agreement within the experimental errors of +25% with those obtained by 
Montgomery. 

Thus, the conclusions, based upon the small values of H^, drawn by 
Hundhausen and Montgomery (1970) are supported by the present measurements, 
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Conclusions 

From these measurements of electron energy fluxes we conclude; 

a) . A detector situated on an interplanetary magnetic field line 
not intersecting the earth’s bow shock measures a net electron energy 
flux away from the solar direction along the field line. The direction 
of this energy flux is independent of the polarity (sign) of the inter- 
planetary field, and very strongly correlated with the instantaneous 
magnetic field direction at all energies up to 10 keV. 

b) . By making comparison between our observations and those of Lin 
and Anderson during a solar electron event, we conclude that the present 
observations extend those of Lin and Anderson to lower energies. Solar 
electrons are apparently emitted with energies as low as 350eV, and contribute 
decisively to during these events. 

c) . When the detector was situated on a line of force of the inter- 
planetary field which had a high probability of intersecting the earth's 
bow shock, it indicated a substantial energy flux at energies >210 eV, 
towards the sun along the field. The energy flux at energies <210 eV 
continued away from the sun. This suggests the possibility of the acceler- 
ation of electrons from energies characteristic of th* plasma up to 
energies of several keV at the shock. The bulk speed of the plasma was 

-i 

~600 km sec at the time of the example given in this paper. 

d) . , the heat flux density, has been measured in four cases by 

calculating up to an energy E^ where the effect of transient electrons 

becomes noticeable. The values obtained, while not of high accuracy, are 

-2 2 

certainly <1x10 ergs/cm /sec, and thus are in agreement with the low 
values measured by Montgomery. The average bulk speed during these 
measurements were approximately 385 km sec~^. 


except on March 19 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


FIGURE CAPTIONS 

Schematic representation of energy flux density <Q^(£)> 
vs. ^ for various physical situations; A), energy flux 
uncorrelated with field direction; B) correlation with 
empty cone about field; C) energy flux filling cone about 
field line. 

Plots of <QqO)>, the differential contributions to the 

energy flux density per unit solid angle, for angular intervals 
o 

of 5 • Note that 3 is the angle between v and B. 

Isometric plots for each detector showing that the largest 
flux is detected by that detector "looking” closest to the 
magnetic field line, and that this condition is preserved 
as the magnetic field direction shifts with respect to the 
triad of directions of the detectors. This is interpreted 
as showing that the energy flux is greatest along the field 
direction away from the sun. 

A plot of <Q^O)> against £ for a period when the interplanetary 

magnetic field intersects the earth's bow shock. The diamond 

shaped points, E u *9.9 keV, refer to the ordinate scale on 

the right, the crosses, E u *210eV, to the scale on the left. 

This figure illustrates a possible difficulty in the 

evaluation of H from measurements of E up to a fixed energy E . 

II II u 

E^ as a function of Log E u for 3 periods during a "delayed” 
electron event (Lin and Anderson, private communication) . 

The crosses (pt. A in the inset) tefer to a time of low solar 
electron flux, the squares (pt. B in the inset) refer to a 
time of high solar electron flux. The event is discussed further 
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Figure 7 


Figure 8 


In the text. Note how extrapolating back to E^O on all 
three curves crosses the axis between 210-340 eV. 

Illustrating how the electron spectra of Lin and Anderson 
join onto those determined by the present experiment. 

Average differential flux values obtained on 8 March and 
11 March 1968 are shown. A correction for the different 
sclld angles of the detectors has been applied, and the 
steepness of the spectrum between 22 keV and 45 keV (power 
law Index -4.4) has been used In Inferring what energy most 
of Lin's particles had. Seventy five percent of the particles 
had energies between 22 and 31 keV. 

The net energy flux E plotted against log E , the uppermost 

II ^ 

energy to which It is calculated. Values of E^X) represent 
flux direction away from the sun, and values <0 represent 
flux toward the sun (away from the earth's bow shock) along 
the magnetic field line. In this example the field almost 
certainly intersects the earth's bow shock. 
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